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Abstract 

The canonical front form Hamiltonian for non-Abelian SU(N) gauge theory in 3+1 
dimensions is mapped non-perturbatively on an effective Hamiltonian which acts only 
in the Fock space of a quark and an antiquark. Emphasis is put on dealing with the 
many-body aspects of gauge field theory, and it is shown explicitly how the higher Fock- 
space amplitudes can be retrieved self-consistently from solutions in the qq-space. The 
approach is based on the novel method of iterated resolvents and on discretized light- 
cone quantization, driven to the continuum limit. It is free of the usual perturbative 
Tamm-Dancoff truncations in particle number and coupling constant, and respects 
all symmetries of the Lagrangian including covariance and gauge invariance. It is 
emphasized that vertex regularization as opposed to Fock-space regularization allows 
an explicit renormalization, which yields linear confinement. The string constant must 
be determined from experiment. 
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1 Introduction 



The Hamiltonian bound-state problem in a field theory is notoriously difficult. Procedures 
like those of Schwinger and Dyson or of Bethe and Salpeter are not easy to cope with in 
practice as reviewed recently in 0]. In the Fifties, the difficulties were considered so enormous 
that the Hamiltonian approach was given up altogether in favor of Feynman's action oriented 
approach. Its modern successors like the lattice gauge theories 0, |[ have arrived at a stage 
of making quantitative predictions In contrast, phenomenological models [H §] are 
closer to experiment and have the different objective to classify the bulk of empirical data. 
They leave little doubt that a heavy meson contains primarily a pair of constituent quarks 
and not an infinity of sea particles as suggested by the perturbative treatment of quantum 
field theory. 

How can one reconcile such models with the need to understand hadronic structure from 
a covariant theory? Particularly, how can one reconcile the constituent quark model (CQM) 
with the quantum field theory of chromodynamics (QCD)? There are several reasons why the 
front form of Hamiltonian dynamics is one of the very few candidates (see fiy, |Tl|, [l^, |T^). 
Particularly the simple vacuum and the simple boost properties [^, [T^, [T^, |T3[ are in contrast 
with the complicated vacuum and the complicated boosts in the conventional Hamiltonian 
theory (see also [|l6l and [17|). Nowadays it is known that even the front form vacuum is not 
simple |p!8| , |T9[, but it is still simpler than in the instant form: The problem can at least be 
formulated mathematically . 



The success of discretized light-cone quantization |]T4|, [T5| particularly in 1+1 dimensions 

H . The two 



has stirred hope that its apparent simplicities carry over to 3+1 dimensions |21 



major obstacles for such an endeavor are the following: (1) The Hamiltonian matrix increases 
exponentially fast with the particle number and the transverse momenta in the Fock states. 
Truncating the Fock space to 2 particles like Tamm and Dancoff |^ invokes perturbation 
theory, violates gauge and Lorentz invariance, and generates non-integrable singularities. Ad 
hoc procedures must be introduced to make things work, both in the original |2^, ^ and the 
light-cone work [Q. (2) The unrestricted integration over the transverse momenta creates 
divergences, which are regularized by some cut-offi The unphysical dependence on the 
regularization parameters has to be removed by renormalization. This is difficult enough in 
a perturbative context, but serious challenges appear when dealing with the non-perturbative 
solutions of a Hamiltonian (see p!7[]). 

In the present work, these problems are addressed and partially solved by adopting 
systematically the point of view that the bound states of QCD are eigenstates of the canonical 
field-theoretical Hamiltonian in the light-cone gauge = 0. The zero modes emphasized 
in previous work are omitted here. In section ^ a method for diagonalizing numerically 
a very large matrix, say with dimension 10^^, by mapping the problem equivalently to the 
diagonalization of a reasonably small matrix, say with dimension 10^, is proposed, using 
the theory of effective interactions E^. The reduction to comparatively small matrices is 



more than a technical issue and requires one to cope with the many-body aspects of a field 
theory. In the present work, the many-body aspects are kept as resolvents which are not 
expanded as perturbative series. They are treated in section ^ by the method of iterated 
resolvents. This novel method requires the development of considerable formal apparatus 
which was given earlier [^. In the continuum limit the matrix equations become integral 
equations, which are similar to those solved numerically already 28|, but here they are 
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derived without assuming that the couphng constant is small and without truncating the 
Fock space. The important issues of vertex regularization and renormalization are discussed 
in sections ^ and ^ lead ultimately to (linear) confinement. The approach is summarized 
and discussed in section ^ 

2 The light-cone Hamiltonian and the Tamm-DancofF 
approach 

In canonical field theory the four commuting components of the energy-momentum four- 
vector are constants of the motion. In the front form of Hamiltonian dynamics P] 



they are denoted by = {P^ , Pj_, P~), see for example |jTT], |T2|. Its spatial components 
P± = {P^jP"^) and P^ do not depend on the interaction and are diagonal operators in 
momentum representation. Their eigenvalues are the sums of the single particle momentum 
components, i.e. 

P^ = T.Pn=^K, and P^ = j:m.- (1) 

Each particle has four-momentum = [p'^ ,p±,p~), sits on its mass-shell p^p^ = and 
therefore has p~ = {rn^ +pf) / p'^ . Each particle state "g" is then characterized by six quantum 
numbers q = {p'^,p±, A; c, /). The first three specify the particle's spatial momentum and A 
its helicity | or |. For QED that's all there is, but for QCD a quark is specified further by 
its color and fiavor, c and /, respectively. Correspondingly, gluons are specified by the five 
quantum numbers q = {p^,Pl, A; a), with a being the glue index. The temporal component 
P~ = 2P+ depends on the interaction and is a very complicated non-diagonal operator. It 
propagates the system in the light-cone time = + x^, i.e. = -P+l^) aiid is the 

proper front-form Hamiltonian 

The contraction of P'^ is the operator of invariant mass-squared, i.e. 

pi^p^ = P+p- - Pi =Hi^^ = H . (2) 
It is a Lorentz scalar and referred to somewhat improperly but conveniently as the 'light- 



cone Hamiltonian i^Lc' shortly H. Its matrix elements (and eigenvalues) somewhat 

unusually carry the dimension of an invariant-mass-squared. In this work one aims at a 
representation in which H is diagonal, 

H\^,) = E,\^^) . (3) 

Since the Hamiltonian does not conserve particle number, the eigenfunctions are com- 
plicated superpositions of particle number projections with contributions like (qql'^i) or 
{qqgl'ifi) ■ ■ ■■ The eigenvalue equation (^ stands thus for an infinite set of coupled integro- 
differential equations which are even difficult to write down explicitly []T2| . 



In order to analyze the structure of the eigenvalue equation, we apply a technical trick, 
referred to as Discretized Light Cone Quantization (DLCQ) 0, |l^ : One imposes periodic 



boundary conditions on the Lagrangian fields. In consequence one maps the infinite set of 
coupled integral equations onto a finite set of coupled matrix equations. The kernel of the 
integral equation (^) becomes a finite matrix, which in principle can be diagonalized by 
numerical methods. 
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Table 1: The Hamiltonian matrix H = T+V for QCD and fixed value of harmonic resolution 
K = 4. — Each element of the matrix below corresponds to a block matrix with potentially 
large dimension, caused by all possible combinations of the single particle momenta subject 
to fixed total value. The diagonal blocks are denoted by T. They are diagonal matrices, 
its elements being the sum of all 'kinetic energies' of the particles in the respective Fock 
states. Due to the selection rules most of the block matrices are zero matrices, marked by 
a dot (■). The blocks marked by V denote the large block matrices which are potentially 
non-zero since the vertex interaction can get active. Their actual matrix elements depend 
on the momenta and helicities of the particles as tabulated elsewhere. 
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Why finite? As a consequence of discretization, the Fock states are denumerable and 
orthonormal. Since has only positive eigenvalues and since each particle has a lowest 
possible value of the number of particles in a Fock state is limited for any fixed value of 
the harmonic resolution K, defined by Eq.(p. One recognizes this best by way of example. 
The lowest possible value K = 1 allows only one gg-pairs (a single gluon cannot be in a 
color singlet). For K = 2, the Fock space contains in addition pairs of gluons, gg-pairs 
plus a gluon, and two gg-pairs, as illustrated in Table |l|. For K = 4, one has at most 8 
particles, namely 4 gg-pairs. In the table all 13 Fock-space sectors possible for K < 4 are 
enumerated by n = 1, . . . , 13. Note that the above classification does not change when K is 
increased further: one just adds more complicated Fock-space sectors. Their number grows 
quadratically with K and has the value Nk = {K + 1){K -|- 2)/2 — 2. 

The single particle momenta can take all possible values subject to the constraints in 
Eq.(|l|). Each of the above Fock space sectors contains infinitely many individual Fock states, 
since the transverse momenta can take either sign. Their number can be regulated by some 
convenient "Fock-space regularization" : A Fock state with n particles is included only if 
the single particle four-momenta satisfy [l^ (pi + P2 + • • - Vuf — (j^i + rn2 + ■ ■ ■ nrinf < 
A|,g. However, it was not realized in the past |T^, that Fock-space regularization is almost 
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irrelevant in the continuum theory. One needs a better alternative, and we propose to work 
with "vertex regularization" . At each vertex, a particle with four-momentum is scattered 
into two particles with respective four-momentum pi and P2- One can regulate the interaction 
by a condition on the off-shell mass Ai'^ = {pi + P2Y — {mi + 7712)^, 

(4) 

in terms of the regulator mass scale A. This innocent step will turn out as a rather efficient 
tool for analyzing the divergences, see below. To get more insight, the relevant momenta are 
parametrized as p^^ = {p^,Pl,P~), Pi = {zp^, zp±_ -|- and P2 = ([1 — z]P^ 1 W ~ Ap^ ~ 

l±,P2)- The off-shell mass becomes then simply 

= (Q + (mi + m2f (z - zf) , (5) 

z[l — z) ^ ' 

and the mass scale governs how much the particles can go off their equilibrium values = 
and 'z = mi/ {mi + 777,2). For A = 0, the phase space is reduced to a point, and consequently 
the interaction is reduced to zero: The Hamiltonian matrix (or the integral equation) is 
diagonal, and has the spectrum of the free theory irrespective of the matrix dimension 
governed by Kps- 

Once the Fock space is defined, one can calculate the matrix elements of the (light-cone) 



Hamiltonian pT| , pT^ . It is written here as a sum of three parts, H = T -\- V + W . The 
kinetic energy T survives the limit of the coupling constant g going to zero. It is diagonal in 
Fock-space representation, its eigenvalue is the free invariant mass squared of the particular 
Fock state. The vertex interaction V is the relativistic interaction per sc. It is linear in g and 
changes the particle number by 1. Matrix elements of V which change the particle number 
by 3 (as in the instant form) vanish kinematically: The vacuum does not fluctuate, see for 
example [jl2|. The instantaneous interactions W are consequences of working in the light- 



cone gauge, are gauge artifacts, and are proportional to g^. They either conserve particle 
number (seagull interactions), or change it by 2 (fork interactions). For the moment we will 
drop the instantaneous interaction W , but we will restore it later. Most of the Hamiltonian 
matrix elements are then zero. Only those where the particle number is changed by 1 are 
potentially non-zero, as displayed in Table |l[ The diagonal blocks are also essentially zero 
matrices, except for the kinetic energies which reside on the diagonal. It is a long standing 
problem of field theory of knowing how such a matrix can have bound states. 

The key problem is that the dimension of the Hamiltonian matrix increases exponentially. 
Suppose, the regularization procedure allows 10 discrete momentum states in each direction, 
i.e. in the one longitudinal and the two transverse directions of A particle has then 
roughly 10^ degrees of freedom. A Fock-space sector with n constituent particles has thus 
10"~^ different Fock states, since they are subject to the constraints, Eq.([I|). Sector 13 alone 
with its 8 particles has the dimension of about 10^^ Fock states. 

The Hamiltonian method applied to gauge theory therefore faces a formidable matrix 
diagonalization problem. Sooner or later, the matrix dimension exceeds imagination, and 
other than in 1-|-1 dimensions one has to develop new tools by matter of principle. One 
would like to reduce the full problem to an effective interaction between a quark and an 
antiquark, 

i/eff|^,5)=^|^,,-). (6) 
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The derivation of such an effective interaction, is the subject of this work. 

Effective interactions are a well known tool in many-body physics In field theory 
the method is known as the Tamm-Dancoff-approach, since it was applied first to Yukawa 
theory for describing the nucleon-nucleon interaction by Tamm pSf and by Dancoff . Let 
us review it in short. 

In a field theory, the Fock-space sectors \i) appear in a most natural way, see Table ||. As 
explained above, the Hamiltonian matrix can be understood as a matrix of block matrices, 
whose rows and columns are enumerated by i = 1,2, ... N in accord with the Fock-space 
sectors like in Table ^ Correspondingly, one can rewrite Eq. as a block matrix equation: 

N 

imj) {m=E foralH = l,2,...,iV. (7) 

One recalls that each sector contains many individual Fock states with different values of x, 
•px_ and A. The rows and columns of the matrix can always be split into two parts. One speaks 
of the P-space P = Z)j=i \ with 1 < n < A^, and of the rest, the Q-space Q = 1 — P. 
Eq. (H) can then be rewritten conveniently 2 block matrix equation 

{P\H\P) {P\^) + {P\H\Q) {Q\^) = E{P\^), (8) 
and {Q\H\P) {P\<i^) + {Q\H\Q) {Q\^) = E {Q\^) . (9) 

Rewrite the second equation as {Q\E — H\Q) {Q\^) = {Q\H\P) (-P|\E'), and observe that 
the quadratic matrix {Q\E — H\Q) could be inverted to express the Q-space wavefunction 
(Ql^) in terms of the P-space wavefunction (P|\E'). But here is the problem: The eigenvalue 
E is unknown at this point. One therefore solves first an other problem: One introduces 
the starting point energy u; as a redundant parameter at disposal, and defines the Q-space 
resolvent or propagator as the inverse of the block matrix {Q\u! — H\Q), thus 

{Q\<^>{u:)) = Gq{u){Q\H\P){P\^>) , with Gq(^) = . (10) 

Inserting this into Eq.@ defines an effective interaction 

H,s{u) = H + H\Q) Gq{uj) {Q\H , (11) 

which by construction acts only in the P-space 

H,s{uj)\P) (P|^fc(c^)) = Ekiu) l^fc(^)) . (12) 

The effective interaction is thus well defined: It is the original block matrix {P\H\P) plus 
a part where the system is scattered virtually into the Q-space, propagating there (Gg) by 
impact of the true interaction, and scattered back into the P-space. Every value of u defines 
a different Hamiltonian and a different spectrum. Varying u one generates a set of energy 
functions Ek{uj). Whenever one finds a solution to the fix point equation 

Ek{uj) = u, (13) 

one has found one of the true eigenvalues and eigenf unctions of H, by construction. 
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One should emphasize that one can find all eigenvalues of the full Hamiltonian H, irre- 
spective of how small one chooses the P-space. Explicit examples for that can be found in 



26]. It looks as if one has mapped a difficult problem, the diagonalization of a huge matrix 
onto a simpler problem, the diagonalization of a much smaller matrix. But this is only true 
in a restricted sense. One has to invert a matrix. The numerical inversion of a matrix takes 
about the same effort as its diagonalization. In addition, one has to vary u and solve the fix 
point equation (p!3|). The numerical work is thus larger than the direct diagonalization. 

The advantage of working with an effective interaction is of an analytical nature because 
the resolvents can be approximated systematically. The two resolvents 

''«''^' = (oi.-r-f/|Q) """^ = {Q|.- no)' 

defined once with and once without the non-diagonal interaction U, are identically related 
by Gq{u!) = Go{uj) + Go{uj) U Gq^u), or by the infinite series of perturbation theory. The 
point is, of course, that the kinetic energy T is a diagonal operator which can be trivially 
inverted to get the unperturbed resolvent Gq^uj). In practice p2|, ^ ^ one has never 



gone beyond the first non-trivial order of perturbation theory for formulating the effective 
interaction, namely He^i^^) = Tgg + W + V\qqg) Go{lj) {qqg\V. This, of course, is highly 
unsatisfactory. 

The identification of the gg-sector with the P-space and the related effective interaction 
appears as the most natural thing to do. But in practice, the Tamm-Dancoff-approach 
(TDA) has two serious defects: (1) The approach is technically useful only if one truncates 
the perturbation series to its very first term. This destroys Lorentz and gauge invariance. 
(2) If one identifies u with the eigenvalue as one should, the effective interaction develops 
a non-integrable singularity, both in the original instant form |2^, and in the front form 



2^, [23]. In order to force the approach to work, one has to replace the eigenvalue cu by the 



kinetic energy in the P-space Tgg [|^, . 



3 The effective interaction and the wavefunction 

In gauge theory the Fock space sectors appear in the most natural way, and in DLCQ the 
number of Fock-space sectors is finite from the outset as shown in Table |1|. The block matrix 
is sparse due to the nature of the underlying Hamiltonian, most of the block matrices are 
zero matrices. One should make use of that! The Tamm-Dancoff approach, particularly the 
step from Eqs.(P) and (P) to Eq.(|l2]), can be interpreted as the reduction of a block matrix 
dimension from 2 to 1. But there is no need for identifying the P-space with the lowest 
sector. Nothing prevents us from choosing the last sector as the Q-space: The same steps as 
above reduce then the block matrix dimension from to — 1. The effective interaction 
acts in the now smaller space. This procedure can be iterated. The price to pay is that 
one has to deal with 'resolvents of resolvents', or with iterated resolvents. Ultimately, one 
arrives at block matrix dimension 1 where the procedure stops: The effective interaction 
in the Fock-space sector with only one quark and one antiquark is defined unambiguously. 
In order to make progress one needs to readjust notation p6[. Suppose, in the course of 
this reduction, one has arrived at block matrix dimension n, with 1 < n < N. Denote the 
corresponding effective interaction Hn{oj). Since one starts from the full Hamiltonian in the 
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Figure 1: The three graphs of the ef- 
fective interaction in the qq-space. The 
lower two graphs correspond to the chain 
U = VG3V, the upper corresponds to Ua = 
VGSVG2VG5V . Propagator boxes are repre- 
sented by vertical dashed lines, with the sub- 
script 'n ' referring to the sector numbers. 



last sector A^, one has to convene that = H. The eigenvalue problem corresponding to 
Eq.([T2|) reads then 

f:(z|if„H|j)(j|M/H) = EM , forz = l,2,...,n. (15) 

Observe that i and j refer here to sector numbers. Now, in analogy to Eq.(|10[), define 

n—l 1 

{nl^ico)) = Gniuj)J2{n\Hnico)\j) (jl^iu)) , with Gniu;) = — ^ , ., , ■ (16) 

p[ {n\uj - Hn{io)\n} 

Note that this definition requires only the inversion of the sector Hamiltonian Hn{uj), but 
not of the total Hamiltonian. The effective interaction in the (n — l)-space becomes then 

Hn^liuj) = Hn{uj) + Hn{uj)Gn{uj)Hn{uj) . (17) 

This holds for every block matrix element {i\Hn_i{uj) To get the corresponding eigenvalue 
equation one substitutes n by n — 1 everywhere in Eq. ([T5D . Everything proceeds like in 
section 0, including the fixed point equation E{uj) = uo. But one has achieved much more: 
Eq. ([T7|) is a recursion relation which holds for all 1 < n < A^! 

The method of iterated resolvents was applied in to the paradigmatic case of a 4 x 4- 
matrix, including explicit numerical examples. Applying the method to the block matrix 
structure of QCD, as displayed in Table |I|, is particularly easy and transparent. By definition, 
the last sector contains only the diagonal kinetic energy, thus if 13 = T13. Its resolvent is 
calculated trivially. Then H12 can be constructed unambiguously, followed by Hn, and so 
on, until one arrives at sector 1. Grouping them in a different order, one finds for the sectors 
with one gg-pair: 

= Ti+ VG3V + VG3VG2VG3V , (18) 
H3 = T3 + VGqV + VGeVG^VGeV + VG4V , (19) 
He = Tg + VGioV + VGioVGgVGioV + VG7V , (20) 
^10 = T^o + [VG,,V + VGi5VGuVGi5V] + VGi,V . (21) 

The quark-gluon content of the respective sectors is added here for an easier identification. 
The square bracket in the last expression indicates how the effective interaction looks for a 
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higher K; for K = 4 the content of the square bracket has to be set to zero. Correspondingly, 
one obtains for the sectors with two gg-pairs: 

H,ggg = H, = T,+ VGrV + VG^VG.VG^V , (22) 
H,g,gg = Hj = Tj+ VG,,V + VGnVG,,VG,,V + VG^V , (23) 
H,^,gga= Hn = Tn + [VG,,V + VG,eVG,,VG,eV] + VG^V , (24) 

and for those with with three gg-pairs 

Hmm = Hs = Ts + VGuV + VGuVGnVGuV , (25) 
H,g,g,gg= Hu = Tu + [VGnV + VGirVGieVGuV] + VGr3V . (26) 

In either case, their structure is different from the one in the pure glue sectors which is 

H2 = T2 + VGsV+ VG5V , (27) 
H5 = T5 + VGeV+ VGgV , (28) 
Hg = Tg + VGioV + [VGuV] . (29) 

Having arrived at this point, one can restore first the limit K — > 00, simply by suppressing 
the square brackets in the above relations. Second, the instantaneous interaction W omitted 
thus far can be restored ex post. Due to the peculiar structure of the gauge theory Hamilto- 
nian displayed in Table |I], the vertex interaction appears only in even pairs, typically in the 
combination VGV. It is plausible that by simply substituting VGV W + VGV restores 
the full Hamiltonian H = T + V + W, as was checked explicitly in [P^]. 

The full effective Hamiltonian H^s = Hqq is given by Eq. (0) and displayed diagrammat- 
ically in Fig. What does it mean physically? The first term is the kinetic energy (Ti) 
in the gg-space. In the second term (VG^V), the vertex interaction V creates a gluon and 
scatters the system from the gg- into the qqg-spa.ce. As indicated in Fig. |l| by the vertical 
line with the subscript '3', the three particles propagate there in '3'-space under impact of 
the full interaction before the gluon is absorbed. The gluon can be absorbed either by the 
antiquark or by the quark, represented by the two respective graphs Ui^i and Ui^2 in the fig- 
ure. The third term {VG3VG2VG3V) , finally, in Eq.(|l^) describes the virtual annihilation 
of the gg-pair into two gluons, and is represented by the graph Ui^q in Fig. |I]. It generates 
an interaction between different quark flavors. As a net result the interaction scatters a 
quark with helicity Xq and four-momentum p = {xP~^, xP± + k±,p~) into a state with A'^ and 

four-momentum p' = {x'P^,x'P± + fc^,j9'~), and correspondingly the antiquark. 

In the continuum limit, the resolvents are replaced by propagators and the eigenvalue 
problem Hes\ip) = M'^lip) becomes again an integral equation but one with a more trans- 
parent structure: 

M^{x, k±; Xq, XqliJb) = 

+ [ dx'£k[ {x, 4; A„ Xq\U + Ua\x', fcj; a;, A;.) {x', k[; X'q, X'ql^b) . (30) 

The effective quark mass m includes the impact of diagram ?7i 2. The domain D restricts 
integration in line with regularization (see below). The effective interaction U includes all 



99 



H, 



999 



H, 



gg 99 



mq + kl ^771^ + kl 



X 



X 



{X, k^] Xq, Xq\lpb) 



9 



fine and hyperfine interactions. The flavor- changing annihilation interaction Ua = t/i,o is 
probably of less importance in a first assault and will be disregarded in the sequel. The 
eigenvalues refer to the invariant mass of a physical state and the corresponding wave- 
function {x, k^; \g, Xglipb) gives the probability amplitudes for finding in that state a flavored 
quark with momentum fraction x, intrinsic transverse momentum k± and helicity A^. Both 
are boost-invariant quantities. Color plays a subordinate role because the Fock states can 
be made color-neutral. 

The eigenfunctions ipb represent the normalized projections of the full eigenfunction 
onto the Fock states \q;q) = 6j(i^|fac). Their knowledge is sufficient to retrieve all desired 
Fock-space components of the total wavefunction |\E'). Let us work that out explicitly by 
asking for the probability amplitude to find a \gg)- or a \qq g) -state in an eigenstate 1^^) of 
the full Hamiltonian. The key is the upwards recursion relation Eq.(|T^). Obviously, one can 
express the higher Fock-space components (n|\E') as functionals of ipqq by a finite series of 
quadratures, i.e. of matrix multiplications or of momentum-space integrations. One needs 
not solve another eigenvalue problem. The first two equations of the recursive set in Eq. (^) 
are 

(2|^) = G'2(2|i/2|1)(1|^), (31) 
and (3|^) = G'3(3|//3|l)(l|^) + G'3(3|//3|2)(2|^). (32) 

The sector Hamiltonians iJ„ have to be substituted from Eqs. ( [TP| ) and ( [27| ) . In taking block 
matrix elements of them, the formal expressions are simplified considerably since many of the 
Hamiltonian blocks in Table ^ are zero. One thus gets simply {2\H2\1) = {2\VG3V\1) and 
therefore (2|^) = G2VG3V{1\^). Substituting this into Eq.(^ gives (3|^) = G3V'(1|^) + 
G3VG2VG3V{1\'^) . These findings can be summarized more readable as 

li^ss) = G,,VGqg,V\^g,), (33) 

and \lPggg) = GqqgV \^ gg) + G gg gV G g gV G gg gV ^^gg) . (34) 

Correspondingly, one can write down the higher Fock-space wavefunctions as functionals of 
only \ipgg). Of course, one has to readjust the overall normalization of the wavefunction, 
and this depends on how many Fock spaces one decides to include. Finally, it should be 
emphasized again that the finite number of terms contrasts strongly to the infinite number 
of terms in perturbative series, because the iterated resolvents re-sum the series to all orders 
in closed form. The above results are all exact. 

Aiming for the higher Fock space amplitudes is not an esoteric problem. One needs them 
for calculating the structure functions, for example. They are the probability amplitudes 
for finding in a hadron a quark with longitudinal momentum fraction x, irrespective of 
its transverse momentum, and to get them, one must sum over all Fock-space amplitudes 
f{x) = J^ni'^nW {x)b{x)\'^n) , scc for example fT^ . 



Let us discuss in some greater detail the structure of the sector Hamiltonians, particularly 
the effective interaction in the qqg-space as given by Eq.(|T9|). The corresponding graphs are 
displayed diagrammatically in Figs. |^ and ^. Those in Fig. |] differ from those in Fig. |l] only 
by an additional gluon. The gluon does not change quantum numbers under impact of the 
interaction and acts like a spectator. Therefore, the graphs in Fig. ^ will be referred to as the 
'spectator interaction' U 3. In the graphs of Fig. |^, the gluons are scattered by the interaction, 
and correspondingly these graphs will be referred to as the 'participant interaction' U3. The 
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Figure 2: The three graphs of the spectator 
interaction in the qq g-space. Note the role 
of the gluon as a spectator. 

Figure 3: Some six graphs of the participant 
interaction in the qq g-space. 



1 




_ 1 _ 


"~ 1 " 


^^^^^ 1 M M 


(33) 4 


(34) 4 

_ 1 _ 




^^^^^^^ \ 


1 




(35) 6 


(36)"" 6 * 




1 


(37) 6 


(38* 6 



analogous separation into spectators and participants can be made in all quark-pair-glue 
sector Hamiltonians: 

i7„ = T„+F„ + ^/„ , for n = 3,6,10,15,... . (35) 

The Hamiltonian is additive in spectators and participants. The spectator interaction t/„ 
therefore can be associated with its own resolvent G 



ni 



G„ = 7^^. (while G.^^— = — ^= (36) 

with an exact relation to the full resolvents Gn, i-e. Gn = Gn + Gn Un Gn- Equivalently, one 
can write it as an infinite series Gn = Gn + GnUn Gn + GnUn GnUn (?„ + ..., as usual. The 
difference to the usual perturbative series in Sec. ^is, that there the 'unperturbed propagator' 
Gq{uj) refers to the system without interactions while here the 'unperturbed propagators' 
Gn contain the interaction in the well defined form of [/„. One therefore deals here with a 
'perturbation theory in medium'. Note that the present series implies also different physics: 
The system is not scattered into other sectors but stays in sector n. This is reflected in the 
operator identity, which is obtained straightforwardly from Eqs.(|35D and (^6]): 



(1 - UnGn){uJ - Hn) = ( 1 - Un ^ ) {uj - H n) = UJ - Hn . (37) 

Taking the inverse gives Gn — GnR^i thus 

Rl = . (38) 

1 — GnUn 

Since Rn and G„ commute, one ends up with Gn = RnGnRn- In all of the above effective 
interactions the square matrices Rn are sandwiched between a quark-pair-glue resolvent G 
and the vertex V , 

VGnV = VRnGnRnV = VGnV, (39) 
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with V = RnV being an effective vertex. The index n needs not to be carried on. One can 
thus rewrite systematically Eqs.([T8|)-(pTD. Working backwards yields 

He = n + VGioV + VGioVGgVGwV, (40) 
Hs = Ts + VGeV + VGeVG^VGeV, (41) 
Hi = Hi = Ti + VG3V + VG3VG2VGSV . (42) 

Instead of being similar, the quark-pair-glue sector Hamiltonians Hn = T„ + ?7„ are now all 
identical. Note that, in the strict sense, they are equal only in the continuum limit. For any 
finite K, there is a last sector which is creating edge effects. The resolvents G2, G5 and Gg 
carry no bar. In the pure glue sectors, the distinction between participants and spectators 
makes no sense. 



4 Analysis of propagators and vertex functions 

Thus far the approach is formally exact. In the sequel, this rigor is given up in favor of 
transparency. For calculating the effective interaction U one needs the resolvents G3. For 
calculating G3, one needs G4 and Gq, and so one. Does this ever come to an end? A careful 
analysis shows that there is a simple relation between them since all Hn{uj) are all bona fide 
Hamiltonians. Suppose we have solved the eigenvalue problem in the gg-space, 

Y^{q;q\H,,ico)\q';q') {q';q'\M^)) = (u) {q; q\M^)) . (43) 

The eigenvalues are enumerated by b = l,2,...cxo. The eigenfunctions (g; g|'?/'f,(u;)) are 
a complete set. Despite working in the continuum limit, we continue to use summation 
symbols for the sake of a more compact notation. Suppose further that an u was found 
which has the same value as the lowest eigenvalue = Mf. The substitution u = 
will hence forward be done without explicitly mentioning. Next, ask for the eigenvalues and 
eigenfunctions in the qqg-space: 

E (g;g;^|:H^3|g';g;^')(g^g';^7VU = ^m(9;^;^I^6,.)- (44) 

It is important to note that one needs not to perform another calculation, becuase the 
result is already known. By construction, the gluon is a free particle which moves relative 
to the meson subject to momentum conservation. The eigenfunction is a product state 
\'4^b,s) = I'^b) ® \ fs)- Parameterizing the gluon's four-momentum as 

q^={yP^,yP± + qL,P;) (45) 

implies that the gluon moves with momentum fraction y and transverse momentum qj_- It 
also implies that the meson moves with momentum fraction 1 — y and transverse momentum 
—q±. Since the gluons are massless, the eigenvalues are therefore 

..2 + ql ql 
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Figure 4: The qq vacuum polariza- Figure 5: The gg vacuum polariza- 
tion graph. tion graph. 

Knowing the eigenvalues and eigenf unctions, one can calculate the exact resolvent. A few 
identical rewritings give 

{q; q; glGsW; 4\ a') = GsW, t, 9) ^; ^k'; 9') - {q; q; 9\^W; q; 9')] , with (47) 

G3{q;q;g) = --^ , and = {y^M' + ql)-^ . (49) 

Since the operator A by definition cannot become a Dirac-5 function {q;q]g\q']q';g'), one 
can drop A, as an approximation. The same type of approximation parametrizes the two 
gluon propagator G2{gi; (72) = {M'^ — M'^b)~^ terms of the glue ball mass M|^. 

Let us digress for a short discussion. One should note first that the above approximation 
can be controlled e posteriori: Once the eigenfunctions in the gg-sector are available numer- 
ically, one can check whether the exact definition of the resolvent Eg. (^71) is peaked like a 
5-function, with a residue predicted by Eg . (^9]) . Second, the above notation is suggestive for 

being related to the single-particle four-momentum transfer along the guark line in graph 
Ui^i in Figure |I[ In fact one gets (p — p'Y = — [|/^(2m)^ + ql], for sufficiently small y. The 
momentum transfer is the same as in Eg.(^) if one replaces the 'current mass' m by the 
'constituent mass' M/2. Apart from that the exact resolvent agrees with the perturbative 
propagator: In the solution, the particles propagate like free particles. Third, and most 
importantly, instead of having resolvents of resolvents, the hierarchy of iterated resolvents is 
broken. Ex post, this justifies the ad-hoc procedures in the Tamm-Dancoff approach, both 
in the original work p^, and in its light-cone adaption p2|, Last but not least, one 
should emphasize that the complicated operator A in Eg. (|47|) can perhaps be omitted when 
dealing with the low energy part of the spectrum. This becomes inadeguate if not false 
at sufficiently high excitations, at sufficiently large u, where one might be in the region of 
'overlapping resonances' |2^. This reguires different technigues, possibly those based on 



random matrix models pO| , |31[| . It is here where the concept of a temperature will possibly 
make its way into the Hamiltonian bound-state problem. 

Next, how can one approximate the vertex functions -R^? With Eg.(|37D one has 



Rn=l + \GJJn + GnU^ + .... (50) 
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For n = 3, this reads R3 = {1 — GsIVRqGqRqV + VRiGiRiV]) 2. Expanding up to the 
first non-trivial term in the expansion gives 



1 + - (Fi + + Tg + r4 + Tg) 



(51) 



For being consistent one must set i?4 = -Rg = 1- Unfortunately, this was not realized in 
the previous work ||26|. The vertex function R^ identifies itself as the conventional vertex 
function. The altogether 5 diagrams of light-cone QCD had been calculated before by Thorn 
and by Perry |Q . Because they are not available in the form we need them in the sequel, 
we have recalculated them and present them elsewhere |3^. In the nomenclature of |Q, the 
mass term Fi cancels almost completely with the two vertex corrections F2 and F3. In the 
sequel we consider explicitly only the two vacuum polarization diagrams F4 and F5, as given 
in Figs. § and Parameterizing the particle momenta as 

Qq = {yzP^, yzP± + zq_i + l_L, p^J , 

f = {y{l - z)P+, yil - z)P^ + (1 - z)q^ - L, p-^) , (52) 
we calculate first the propagators in the qq qq and qq gg sectors, 

g+ zil - z) 



mm 



G. 



mgg 



Gil 



Gp, 



P+ 2(1 - z)Q2 + il + m) 

Qg Z{1 - Z) 

'P+ z{l~z)Q^ + ll ■ 



and 



(53) 
(54) 



They depend on Q^, the momentum transfer along the quark line. The vacuum polarization 
diagrams are then calculable according to the rules of light-cone perturbation theory |jl2[ 
and yield 



^ a 

aric 



dzdh 9o 



dz dl I Qr 



1 - 2z{l 



2mi 



z - 



Z5_ + mj 
2 - z(l 



z{l-z)Q^ 



1 



z{l-z) 



ll + z{l-z)Q'' 



(55) 
(56) 



The dependence on is caused by the propagators. The term with 1/-2(1 — z) in F5 cancels 

. These authors have however not kept 



with Fi + F2 + F3, like in the previous work [32, [ 
track of the Q^-dependence. Dropping this term replaces F5 with F5. The cut-off functions 
are defined by vertex regularization in Eq.(0) and (13) 



QizAi 



1, for < /I < lliz) and < z < 1 - e. 



0, otherwise. 

The limiting transverse momentum l\{z) describes a semi-circle in the appropriate units 



(57) 



lliz) 



with h 



(A^ + (mi + 

+ 2mi(mi + m2 
A^ + (mi + m2)^ 



- - \ z- - 



and 



A^ 



A^ + 4mim2 



A^ + (mi + m2)^ A^ + (mi + m2)^ 



(5J 
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The semi-circle is limited by its intersections at zi = ei and Z2 = I — eu, with 

e; = ^-y^ and e„ = 1 - . (59) 

For sufficiently large A and equal masses mi = m2 = m, the limits become ei = eu = (m/A)^. 
For the purpose of regularization, the gluons are endorsed with a small mass rrig, and thus 
vertex regularization regulates all divergences on the light cone. One gets straightforwardly 



r. = -i:j^ /d.-(i-2.(i-.)-y) m . (60) 

This is approximated by replacing the term m^/2;(l — z) by its maximum at z = ^ that is 



by Attij. The same type of approximation is done for Fs, yielding finally 



F4 = y — — T- - 1 n 1 + . , and F5 = n 1 + . 61 

In either case, the terms depending on e/ = {rrif/Ay or eg = (m^/A)^ have been dropped 
since they are small. The effective coupling constant = aR^ = a{l + F4 + F5), 



, lla^n. , / A^ \ ^ /Sm^, \ , / A^ \ , , 

=.W^)^a + ^l.>(l + 5^)+g^i:(^-l)ln(l + j^^). (62) 



depends therefore on the momentum transfer. Finally, the effective quark and gluon masses 
become 

2 2 2 '^r ~ 1 1 2 2 2 1 / A^ \ 

mr = m. + mr m — ^, and m„ = m„ > rrif m 1 H (63) 

^ ^ ^TT 2rac m2 ' 3 47r ^ V 4m^y ^ ' 



They are obtained straightforwardly by light-cone perturbation theory |]T2[ combined with 
vertex regularization. The diagram for the effective quark mass = mj(l + Fg) is given 
in Fig. ^. There is a corresponding diagram for the effective gluon mass m^. 

5 Renormalization group analysis and confinement 

Collecting terms and restricting to nc = 3 colors, the effective interaction as defined in 
Eq. (|30|) becomes without the flavor changing interaction Ua 



E (x, kx_] Xg, A, 



+ 



X 1 — X 



{X, k±; \g, Xg 



A' 



Y: [ dx'd% ^ (A„ Xg\S{Q)K (^'' ^i' ■ (64) 



One recalls that the four-momentum transfer along the quark (or anti-quark) line is = 
— {kg — k'^y = —{kg — k'^y ~ {nig + mqY{x — x'Y + (A;j_ — fc[)^. As we are supposed to 

15 




Figure 6: The quark mass diagram. 



m„ = 


6 MeV 


rrid = 


10 MeV 


rris = 


200 MeV 


rric = 


1.40 GeV 


rrib = 


4.70 GeV 


m = 


180 GeV 


a = 


0.2 


A = 


100 MeV 




13 MeV 



Figure 7: Fixing the bare quantities. 



do, we restore the instantaneous interaction W in the gg-space when deriving the above 
result, thus Ugg = Wgg + V GqqgV . Both the instantaneous interaction W and the vertex 
interaction V GgggV have a non-integrable singularity, c/]T2|, which cancels each other such 
that only the integrable Coulomb singularity remains. The latter manifests itself in Eq. (|6^) 
in the "Coulomb denominator" Q~^. The spinor factor S{Q) represents the familiar current- 
current coupling which describes all fine and hyperfine interactions 



{Xg,Xg\SiQ)\X'X') 



u{kg, Xg) Y u{k'g, X'g) u{kg, Xg) 7^ u ( fc^ , 



'xx' 



sjil-x)il-x') 



(65) 



One can simplify the expression technically by suppressing in the latter all momentum de- 
pendence. Evaluating the spinors at equilibrium, i.e. at vanishing transverse momenta and 
at momentum fractions Xg = nig/ {nig + rriq) and Xg = nig/ {nig + rriq), gives |1^ 



Uikg, Xg) U{k'^, X'q)] U{kg, Xg) 7^ U{k'g, X'gj 



XqXq 



XqXq 



equilibrium 



4(mg + nig)'^ Sx„x'^Sx,,x'^ • (66) 



Omitting the spin-dependence of the interaction implies to restrict oneselve to the central 
part of the interaction. In this approximation Eg . (|64D becomes 



E (x, ki 



+kl ^ m? + kl 



X 



1 — X 



{x,k_i 



A{mg + mqf 



dx d k I ~ (^x ^ k I 

D 



(67) 

It is this equation, finally, which we are going to analyze by the renormalization group. 

In order to do calculations, one has to specify 10 "input parameters". As specified in 
Fig. ^, these are the six flavor masses ruf the coupling constant a = g'^/4:7i, and the two 
regularization parameters A and rrig. To be concrete, these parameters are specified in Fig. 0, 
and subject to future change. The eigenvalues of the effective equation ([67| ) depend on all 
of them, i.e. Ei{mf, a; A, rrig). While their dependence on the masses and coupling constant 
is not unexpected, they should not depend on A, since that is an arbitrary and unphysical 
parameter. One therefore must require that the eigenvalues satisfy the renormalization group 
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equations 

— Ei{mf,a;A,mg) = 0. (68) 

Since the Lagrangian 'bare' parameters themselves are arbitrary, one can try to let them be 
functions of A in order to satisfy Eq.(|68|): ruf = m/(A), a = a{A), and irig = mg{A). The 
numbers in Fig. |^ can be viewed as such functions at the particular value A = 100 MeV. If 
one finds one single set of functions which makes all eigenvalues Ei independent of A, one 
has solved the problem. It is not really necessary that A — oo. 

In the present work, the eigenvalues of the full light-cone Hamiltonian H^c in Eq. (|^) are 
conceptually identical with the eigenvalues of the effective Hamiltonian, Eq.(^). One must 
therefore analyze the latter. The dependence on the cut-off resides in three well separated 
locations: (1) in the domain of integration D, (2) in the effective coupling constant as{Q, A), 
and (3) in the effective mass parameters rnq{A) and mq{A). As an advantage of vertex 
regularization, one can decouple the scales in D and in the functions as and friq. In the 
sequel we disregard the dependence on D. 

If Ug and frig do not dependent on A, the eigenvalues cannot depend on A, either. It 
therefore suffices to require 

Sas = ^Sa + ^6A' + J2p^6m}^0, and (69) 
da oA^ ^ omj ■' 

^ da dA^ dm) ■> ^ ' 

The variations of Hg have to be performed of course at a fixed value of . Solving the 
corresponding differential equations for a (A) and mj(A) yields the required functions. There 
are two regimes which are particularly interesting. One is the asymptotic regime 



mJ,mJ<Q2^A2, (71) 

in which Eqs. (|6lD and (^) become 

2 1 . , 33 - 2nf 

as = a + a oq In — , with oq = — — — - . (72) 

Solving Eq.(|69D in this regime yields asymptotic freedom in the front form, as shown first by 
Thorn and later by Perry et al. |]53[ . The asymptotic regime is particularly important for 



high energy scattering. But in a bound state calculation, the very large momentum transfers 
are less important than the very small ones, because of the Coulomb singularity ~ 0, see 
Eqs.(^^ and (0). It provides the binding. In the sequel we therefore emphasize the hound 
state regime near ~ 0, 

< < A^ < (73) 

To our recollection, such a regime has not been analyzed by the renormalization group thus 
far. We shall differentiate between the two light quarks with mass mi (u,d) and the four heavy 
quarks with mass m^ (s,c,b,t). For convenience we collect some mass parameters consistent 
with ||3^ in Fig. ^ In the present context, 'light' and 'heavy' is defined operationally by 



mf , mi , , 

and ^»1. (74) 
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mass nig = 0. This in turn defines nig 



Consider first the effective gluon mass. In a gauge theory the vector bosons should have zero 
™?-£i:"^n(l + 3^) + £i:"?ln(l + ^). (75) 



according to Eq. ( ]63| ) . Expanding the logarithm for the former gives 



m 



^ = -A^ f 1 + E ?f In fl + ^) U -A^ . (76) 
^ 47r V Y 4mfy; 47r ^ ' 



With the parameter values of Fig. the numbers in the parenthesis go like (1 + 0.017). It 
is thus save to neglect the impact of the light quark masses. This yields rrig ~ 13 MeV as 
quoted in the figure and removes the regulating gluon mass as a free parameter from the 
theory. In the vicinity oi ~ 0, the effective coupling constant goes like 

„ 2 , / A2 \ 33^2 , / A2 



Gn^^^'V^ Am}) Stt^ + Am} ' ^^^'^ 



The last term is due to expanding the coefficient of up to second order. Treating differently 
the heavy and the light quarks again and substituting the gluon mass gives 

Near ~ the renormalization group equation for the effective coupling constant ( |69D 
gives 



6a 



aA' 



a^A 

—— (5A = , thus 5{aA) = . (79) 



Quite unexpectedly, (a A) turns out as a renormalization group invariant. The coupling 
constant "runs" like 

a{A) = ^ . (80) 

It is admitted to call this "asymptotic freedom in the bound state regime". The subscript 
()o denotes the renormalization point. A particular "renormalization point" is the set of 
numbers in Fig. [7|. Finally, with the same approximations, the effective quark mass becomes 
according to Eq.(p3|) 

_2 nf. 4a, An^ 

m f = m /. 1 H m — 

•' V 37r a / 

Renormalizing the quark mass according to Eq.(|70D and introducing t(A) = In {An /a) as an 
integration variable, thus t{A) = In (47rA/[aoAo]), gives 

ln!!^= fdt . (82) 



mf(0) J t + ^ exp t 
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One notes that the mass ratios mf/mfi = Wif/mfi are renormahzation group invariants, as 
they should. 

In the bound state regime, the kernel of the integral equation (^) becomes thus 

^-Q^ + Q^ + ---- ^^^^ 

Its Fourier transforms behave like (3/r + /t^r, i.e. like the superposition of a Coulomb and 
a linear potential governed by a "string constant k", see for example [^. The theory 
"confines". The two coefficients 

K = ^, and /5 = ao + --illn 1 + - l + -ln^ , 84 

V27r 127r V ao/ 87r ^ V 3 2mi/ 

are related to each other by the renormahzation group. The string constant becomes a 
parameter of the theory which must be determined by experiment. The value k ~ 10 MeV 
consistent with Fig. |^ is probably not very realistic. 

Although the present approach has yet to be applied to hadrons, the calculational scheme 
is rather different from the one taken by Wilson and collaborators |1^, see also [36, 37] and 
mostly because the present approach emphasizes so strongly the region around 

~ 0. A detailed comparison should be given separately when dealing with renormalizing 
the scales connected to the integration domain D in Eqs. ([6^ ) and (jG^). It is more than 
likely that the similarity transform of Wegner and collaborators [^, ^ will be there 
more than useful. 

In order to do calculations, one must give 8 parameters: the bare coupling constant, the 
renormahzation scale and the six bare flavor masses (a,A,m/). As always in the sequel, we 
understand these parameters as given at the renormahzation point and drop the subscript 
Oq. It is not relevant in this context whether one gives the bare or the dressed masses. 



since they are related by the coupling constant in a simple way, Eq.(|81|). These are seven 



dimensionful and one dimensionless parameter. We like to think of them in an other way. 
The relation between A and k is linear, see Eq . (|8^) . Because of that one can introduce 
{a,K,mf) as independent variables and drop any explicit reference to A as the parameter 
of regularization. Rather the string constant k, sets the scale of the problem. The rest are 
dimensionless parameters, the coupling constant or mass ratios, which are introduced to 
cope with the plenitude of experimental numbers. It looks as if the effective theory has one 
parameter more than the Lagrangian, but this is not quite true. In a Lagrangian theory, the 
absolute values of the spectral eigenvalues Ei are physically meaningless. One can multiply 
the Lagrangian and thus the Hamiltonian with an arbitrary number without changing the 
physics. In a way, renormahzation theory fixes that number by explicitly introducing a scale. 

Finally, one should emphasize that one is not completely free in choosing the scale: The 
above formalism holds only for at least one 'heavy' quark, with a mass much larger than k. 
Although there is a long way until one reaches the 180 GeV of the top, one may ask what 
happens to the above when one transgresses this limit. If one repeats the above analysis for 

^ but in the regime -C -C A^ the formalism changes. The situation is similar to 
the asymptotic regime: The coefficient of Q"'^ is roughly m'ja'^, thus unstable with respect to 
the renormahzation group. As a consequence it fades away. The result of a renormahzation 
group analysis depends on the scale at which it is performed. In conclusion, we are afraid to 
pretend that the spectrum of toponium will be that one of a de-confining system, with an 
ionization threshold like positronium. 
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6 Summary and discussion 



The present work has four important ingredients: (1) The method of iterated resolvents 
is apphed to a gauge theory Hamiltonian in the front form; (2) In the solution the full 
quark-pair-gluon propagators can be well approximated by free propagators with the quarks 
carrying the constituent mass m — M/2; (3) Vertex as opposed to Fock space regulariza- 
tion is more efficient to regulate all divergences in the light-cone formahsm; (4) After a 
renormalization group analysis, the spectrum of the Hamiltonian is independent of the two 
regularization parameters and in the confining regime a new mass scale appears which must 
by determined by experiment. — After the dust has settled, one seems to get a glimpse of 
an idea how confinement could work in practice. 

One should emphasize that in this non-perturbative approach there is no Fock space 
truncation. One calculates Lorentz scalars, the eigenvalues of the invariant mass-squared, in 
a frame independent way. Even the kernel of the interaction a(Q)/Q'^ is expressed in terms 
of an invariant, — — (p' — p)"^, which is the momentum transfer of the quark. 

The present work adopts systematically the point of view that all properties of a La- 
grangian gauge field theory are contained in the canonical front-form Hamiltonian, which 
is calculated in the light-cone gauge without the zero modes. Periodic boundary conditions 
allow one to construct explicitly a Hamiltonian matrix. After regularization, one is con- 
fronted with the diagonalization of a large but finite matrix. The binding of the particles 
in a field theory is provided by the virtual scattering into the higher Fock spaces and the 
Tamm-Dancoff approach accounts for that, however it accounts for that only in lowest or- 
der. For large values of the coupling one loses control. These difficulties are overcome in 
the the present work after re-analyzing the theory of effective interactions. By inverting the 
hierarchy, it was possible to develop a new formalism, the method of iterated resolvents, 
and using this method, one can map the large Hamiltonian matrix onto a matrix which is 
smaller. Finally, the periodic boundary conditions are relaxed by going over to the contin- 
uum limit. Having found the eigenfunctions of the effective interaction, one can retrieve all 
other many-body amplitudes in a self-consistent way. 

We saw that the effective interaction has two parts: The flavor conserving interaction 
U and flavor changing interaction Their diagrammatic representations look like second 
order diagrams of perturbation theory, but represent a re-summation of perturbative graphs 
to all orders. Particularly, U bears great similarity with a one-gluon-exchange interaction. 
In analyzing this effective interaction one uses the apparent self-similarity in a gauge theory. 
Ultimately, this aspect allows the 'breaking of the hierarchy': In the solution, the particles 
propagate like free particles. All many-body aspects reside in the vertex coupling functions, 
in which the well known divergences of fleld theory appear typically as kl/x. This can 
become large cither for large transverse momentum k±or for small longitudinal momentum 
fraction x. In order to remove these divergences, one needs two regulator scales. Regulating 
the theory on the level of the fundamental interaction, we introduce vertex regularization 
associated with a mass scale A and a small kinematical gluon mass rrig. Insisting that the 
gauge bosons have vanishing dynamical mass, relates nig to A and removes it from the 
formalism. 

The appearance of a new scale k, provides the room for an interesting speculation. One 
knows empirically that the invariant mass square of the light mesons like the pion is linear 
in the quark masses. For dimensional reasons one needs a parameter of dimension mass 
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to account for that. One believes that this parameter is the gg- vacuum condensate (^^), 
i.e. Ml ~ [niq + mq)(\I/\I/). The vacuum condensate is extremely difficult to calculate and 
usually is taken from experiment. With the appearance of the additional scale which can 
play that role, one can have M^^ ~ (m^ + mq)K. Also k has to be fitted to experiment, but 
its role within the theory is well defined. 
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